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Orientational  Pair  Correlation,  Polar  Order  Parameter  and 
Temporal  Stability  in  the  Second  Harmonic 
Generation  of  Amorphous  Polymers 


Abstract 


Measurements  of  the  second  harmonic  generation  (SHG)  signal  were  made  on  the 
DAN/PMMA  guest/host  system  as  a  function  of  the  nonlinear  optical  (NLO)  chromophore 
concentration  and  the  poling  field  strength.  The  temporal  stability  of  the  SHG  signal,  after  the 
poling  field  was  removed,  was  also  investigated.  The  polar  orientational  order  parameters,  L33 
and  L33  (see  text  for  definition),  were  deduced  from  the  second  order  nonlinear  susceptibility. 
Both  of  the  polar  orientational  order  parameters  and  the  relaxation  time  of  the  nonlinear 
susceptibility  were  found  to  increase  with  increasing  NLO  chromophore  concentration.  The 
concentration  dependance  was  interpreted  as  due  to  the  anisotropic  intermolecular  interaction, 
particularly  the  dipole-dipole  interaction.  The  anisotropic  intermolecular  interaction  introduces 
the  orientational  pair  correlation  function  which  enhances  the  polar  order  as  well  as  the  temporal 
stability  of  the  NLO  effect. 


I.  Introduction 


The  second  order  nonlinear  optical  effect  in  the  amorphous  polymer  lies  in  the  existence 
of  the  polar  orientation  order  induced  by  the  externally  applied  electric  field.  The  external 
electric  field  aligns  molecular  chromophores  which  possess  a  large  hyperpolarizability.  This  type 
of  polymeric  material,  while  combining  the  processability  and  excellent  transparency 
characteristics  of  an  amorphous  glassy  polymer,  represents  a  new  class  of  optical  materials  that 
show  a  great  promise  for  the  development  of  electro-optical  devices  and  second  harmonic 
generation  applications.^’^ 

To  maximize  the  nonlinear  optical  (NLO)  effect  it  is  important  to  understand  the  behavior 
of  the  polar  order  parameter  (POP)  in  association  with  the  effect  of  the  electric  field  poling 
process.  Several  methods  have  been  employed  to  develop  materials  with  large  second  order 
optical  susceptibilities  and  enhanced  stability.  The  incorporation  of  the  NLO  chromophore 
molecules,  by  covalently  bonding  them  to  polymer  chains,  greatly  restricts  the  reorientational 
mobility  of  the  chromophores  and  significantly  improves  the  NLO  stability.  The  covalent 
bonding  of  the  NLO  chromophores  to  polymers  can  be  incorporated  at  the  side  chain,^  main 
chain^  or  through  crosslinkage.^  For  the  chromophore-doped  polymer  guest/host  system, 
relaxation  is  rapid  and  second  harmonic  generation  efficiency  is  generally  short-lived.^ 
However,  due  to  its  ease  of  sample  preparation  with  well  controlled  chromophore  concentration, 
the  chromophore  doped  material  provides  a  well  defined  system  to  delineate  various  physical 
mechanisms  that  affect  the  alignment  and  stability  of  the  NLO  chromophore/polymer  system. 

In  this  paper  we  are  concerned  with  the  concentration  dependence  of  the  POP  and  the 
relaxation  of  the  POP.  The  concentration  dependence  of  the  POP  and  the  relaxation  time  of  the 
POP  is  interpreted  as  being  due  to  the  pair  orientational  correlation  arising  from  the  anisotropic 
intermolecular  interaction  of  the  NLO  chromophores. 
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II.  Experimental 


Appropriate  amountsof4-(N,N*-dimethylamino-3-acetamidonitrobenzene)DANand  (Poly 
(methyl)  Methacrylate)  PMMA  were  dissolved  in  chloroform  to  form  solutions  of  different 
chromophore  concentrations.  The  volume  of  chloroform  in  each  solution  was  adjusted  to  give 
a  desired  viscosity  suitable  for  spin  coating.  The  solutions  were  filtered  to  remove  undissolved 
particulates.  Films  were  prepared  by  spin  coating  the  polymer  solution  on  soda  lime  glass  slides, 
which  were  precoated  with  300  A  Si02  and  250  A  ITO  (indium  tin  oxide)  films.  The  NLO 
polymer/ITO  sample  assembly  was  placed  in  a  vacuum  oven  first  at  40“C  for  over  24  hours  and 
then  baked  at  the  temperature  slightly  above  the  glass  transition  temperature.  The  absence  of 
solvent  was  checked  with  infrared  spectrometry  after  the  baking  process,  no  solvent  IR 
absorption  band  can  be  detected.  After  certifying  the  chromophore  concentration,  another  ITO 
glass  slide  was  then  placed  on  top  of  the  polymer/ITO  glass  slide  to  form  a  sandwich 
configuration  for  electrode  poling.  The  glass  transition  temperature  (Tg)  of  the  sample  was 
determined  by  using  a  DSC  (Differential  Scanning  Calorimeter,  Perkin  Elmer  Delta  series).  The 
refractive  index,  and  the  thickness  of  the  sample,  were  determined  by  a  prism  coupler  (Matricon) 
modified  for  the  multiple  wavelength  operation. 

The  second  harmonic  generation  (SHG)  was  carried  out  using  an  apparatus  similar  to  that 
reported  by  Guan  and  Wang^.  It  consists  of  an  NdrYAG  laser  (Spectra-Physics  GCR-11,  X  = 
1.06  pm,  Q  switched  at  10  Hz),  a  computer  controlled  goniometer  stage,  a  series  of  short-pass 
and  long-pass  filters,  a  polarizer,  a  half  wave  plate,  a  photomultiplier  tube,  and  a  boxcar 
integrator  which  was  interfaced  to  a  PC.  The  sample  assembly  mounted  on  the  goniometer  stage 
was  placed  in  an  oven  which  is  temperature  controlled  to  the  accuracy  of  ±  0.5 °C.  The  electric 
field  poling  was  carried  out  inside  the  temperature  controlled  oven.  To  obtain  the  signal  decay 
measurement  the  electric  field  was  switched  off  and  the  electrodes  were  grounded,  after  the 
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sample  signal  had  reached  a  stabilized  steady  state  in  a  constant  temperature  environment. 

The  Maker  fringe  of  a  Y-cut  single  crystal  quartz  plate  (du  =  0.5  pm/V)  was  used  as 
the  reference  to  determine  the  second  order  susceptibility  of  the  sample. 

m.  Results  and  Discussion 

Shown  in  Fig.  1  are  the  temporal  characteristics  of  the  boxcar  SHG  output  signals  for  a 
DAN/PMMA  film  (DAN  is  10  wt%)  at  84°C,  poled  with  various  levels  of  the  electric  field 
amplitude.  The  sample  assembly  that  was  used  to  obtain  signals  given  in  Fig.  1  was  oriented 
with  the  angle  of  incidence  <f>  =  50°.  The  incident  laser  beam  was  at  the  p-polarization.  At  this 
orientation  the  SHG  signal  was  at  its  maximum  value  with  respect  to  the  SHG  signal  versus  4> 
plot. 

Several  observations  can  be  made  in  Fig.  1:  1)  The  SHG  signal  increases  steadily  with 
increasing  poling  field  strength;  2)  The  leading  edge  of  the  SHG  signal  is  sharp;  3)  considerable 
signal  fluctuations  are  noted  at  the  maximum  amplitude;  4)  After  the  poling  field  is  switched-off 
the  signal  falls  off  rapidly,  and  then  is  followed  by  a  slower  decay.  We  have  investigated  these 
features  and  now  describe  the  result  in  more  detail. 

The  SHG  signal  (Ishg)  investigated  as  a  function  of  the  poling  field  strength  at 

various  NLO  chromophore  concentrations.  Shown  in  Fig.  2  is  the  plot  of  the  SHG  signal  (for 
10  wt%  DNA  in  PMMA  at  84 °C)  as  a  function  of  Ep^,  the  square  of  the  strength  of  the  electric 
field.  To  determine  the  value  of  Ep,  we  used  the  film  thickness  data,  determined  from  the 
waveguiding  experiment  using  the  Matricon,  which  simultaneously  provides  the  refractive  index 
and  the  film  thickness  data  (film  thickness  is  3.89  (xm  for  the  10  wt%  film).  In  the  0.6  to  1.3 
X  10^  volt/cm  range,  appears  to  vary  linearly  with  the  poling  field  as  Ep^. 

For  an  isotropic  film  subject  to  a  poling  field  in  the  direction  perpendicular  to  the  film 
surface,  the  transmitted  second  harmonic  intensity  12^,  ^shg)  excited  with  the  incident 

Page  5 


fundamental  beam  in  the  p  polarization  is  given  by 


I,.  ^  ^  x2(^) 


T2J4>) 


where  c  is  the  velocity  of  light  in  vacuum;  Iq  is  the  transmission  coefficient  of  the  second 
harmonic  light  through  the  substrate.  is  the  intensity  of  the  fundamental  beam  inside  the 
medium.  T2.  {4>)  is  the  transmission  factor,  which  depends  on  the  polarization  of  the  incident 
beam,  ^  the  angle  of  incidence,  as  well  as  the  refractive  index  of  the  material  at  the  fundamental 
(nj  and  doubled  (02^^)  frequencies.  We  have  modified  the  calculation  of  T2^  {<i>)  given  by 
Jerphagnon  and  Kurtz*  for  the  sample  geometry  appropriate  for  the  contact  electrode  poling  and 
corona  poling  techniques.^  The  12^  i<i>)  expression  appropriate  for  the  present  experiment  is 

T  (4))  =  ^°2<.<^os(t)3j(n^cos<{)  ^  %cos^J)  (ja^^cos^J  +n^cos^^J)  ^2) 

(n^^cosil)  +nocos{l)2^')^  (n2„cos<l>  +noCos(|)2<,') 


where  n^  is  the  refractive  index  of  the  front  and  back  substrates;  (f>^^ '  and  4>2w>'  are  the  angles 
of  refraction  of  the  fundamental  and  the  second  harmonic  beams  inside  the  medium,  respectively. 
They  are  related  to  the  incident  angle  <i>  by  Snell’s  law:  sin  (f>^  =  riQsin  4>/n^  and  sin  4>2w'  ~ 
ngsin  4)1^2^.  The  quantity  1^.  is  the  correlation  length  given  by  I,.  =  X/4  ]  n^cos^) J-n2^cos4>2«*  |  ■ 
Finally,  the  effective  optical  susceptibility  (xeff)  for  the  optical  field  in  the  plane  of  incidence  (in 
the  p-polarization)  is  given  by 

=  lx3i^^^sin<^^cos^„'cos^2«'  +  (X33®sin^«i^'  +  x^^^iCos^<^^')sm^2«'|  (3) 

In  Eqs.  (2)  and  (3),  we  follow  ref.  8  by  designating  the  quantity  inside  the  NLO  medium  by  a 
prime. 

The  two  SHG  susceptibility  elements  Xa/^^  and  xaa^^^  are  related  to  the  0  polar 
orientational  order  parameter  1.31  =  Vi  (<cos  >  -  <cos^  >)  and  L33  =  <cos^  0j  >  by^ 
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and 


X31®  -  p  0  flu 


(4) 


X33®  =  P  ^  fix  tj  L33  (5) 

where  p  is  the  number  density  of  the  chromophores  and  /3  is  dominant  component  of  the  second 
order  hyperpoiarizability  tensor  along  the  molecular  axis.  In  the  Lorentz  model  the  local  field 
factors  f2^  and  4  are  given  by  fiu  =  +  2)/3  and  f^  =  (n^^  +  2)/3. 

In  the  POP  expressions  the  subscript  i,  associated  with  d,  refers  to  the  representative 
NLO  chromophore  1. 

For  a  system  of  non-interacting  NLO  chromophores,  the  two  POP’s  are  given  by 

L33  =  c£^(a)=(4  +  l)cotha-(4  +  -)  (6) 

a 

and 

L31  =  i [se^Ca)  -<£3(3)]  =:  (1  +  i.)  - Icotha  (7) 

2  a^  a  a’ 

where  j  and  are  the  Langevin  function  of  order  1  and  3,  respectively.  The  quantity  a  is 
equal  to  Ep/kT,  fg  is  the  local  field  factor  introduced  to  account  for  the  presence  of  the 
polymer  host,  p.  is  the  ground  state  dipole  moment. 

In  the  independent  dipole  reorientation  model,  the  orientation  of  the  dipoles  is  completely 
dictated  by  the  external  field,  and  the  polar  order  parameters  Lgj  and  L33  are  practically 
independent  of  the  concentration  of  the  NLO  chromophores,  due  to  the  fact  that  in  this 
approximation  the  potential  energy  of  the  interaction  of  the  NLO  molecules  does  not  play  a  role 
in  affecting  the  values  of  POP. 

In  the  case  of  the  poling  parameter  a  <  <  1,  Lgj  and  L33  reduce  to: 
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Thus,  under  this  simplified  case  the  NLO  susceptibility  X33^‘^  (or  X3i*^‘’)  is  proportional  to  the 
poling  field  Ep;  hence,  the  SHG  intensity  is  proportional  to  Ep^,  consistent  with  the  result  shown 
in  Fig.  2.  The  ratio  of  X33^"^/X3/^^  is  equal  to  3. 

We  have  determined  the  NLO  susceptibilities  X3i^^^  and  X33^“^  by  fitting  to  Eq.  (1)  the 
measured  intensities  versus  4>  data,  with  the  incident  fundamental  beam  in  the  p  polarization, 
assuming  X33^^^/X3i^"^  =  3. 

An  independent  check  of  the  fitting  procedure  was  also  carried  out  without  making  the 
X33^~Vx3/^^  =  3  assumption,  but  using  the  measured  SHG  intensities  obtained  by  the 
fundamental  beam  in  both  the  s  and  p  polarizations,  a  method  similar  to  that  described  in  ref. 
7.  Within  experimental  uncertainty  the  susceptibility  results  are  in  agreement  with  that  obtained 
by  the  X33^“Vx33^"^  =  3  assumption. 

Having  obtained  the  X33^^^  arid  Xsi^^^  values,  we  can  calculate  the  POP  L33  and  L31,  using 
Eqs.  (4)  and  (5).  We  have  found  that  both  L33  and  L3j  depend  on  the  NLO  chromophore 
concentration.  This  result  is  consistent  with  that  found  in  the  MNA/PMMA  (MNA,  2-methyl  4 
nitroaniline)^,  but  is  not  in  agreement  with  that  predicted  by  the  independent  dipole  orientation 
model  described  above. 


In  Fig.  3,  we  have  plotted  L33  as  a  function  of  the  DAN  concentration  (in  number 
density).  The  /3-value  that  was  needed  for  calculating  L33  (c.f.  Eq.  (4))  was  obtained  from  the 
paper  by  Eaton^®.  One  notes  that  over  the  concentration  range  of  1  x  10^®  to  3.7  x  10^^ 
chromophores/cm^,  L33  increases  from  0.28  to  0.38. 

As  recently  shown  by  Wang,^'  the  concentration  dependence  of  POP  is  due  to  the 
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orientational  pair  correlation  that  arises  from  the  angular  dependence,  or  anisotropic 
intermolecular  potential  between  NLO  chromophores.  In  ref.  11,  it  is  shown  that  while  any  type 
of  anisotropic  intermolecular  potential  (whether  it  be  short  or  long  range  interaction)  can  affect 
the  POP  induced  by  the  external  electric  field.  The  dipole-dipole  interaction  makes  the  most 
important  contribution  because  of  its  long  range  interaction. 

In  the  case  of  the  weak  field  poling  condition  (a  <  <  1)  Wang  shows  that 

L33  =  a  (1  +  p  GJ/5  (10) 

and 

L31  =  a(l  +  pG^)/15  (11) 

where  G^  is  the  cluster  integral  associated  with  the  solution  of  the  molecular  pair  correlation 
function.  In  the  case  of  the  dipole-dipole  interaction,  the  cluster  integral  can  be  related  to  the 
Kirkwood  g-factor,  gj^,  that  arises  from  the  dielectric  theory  of  the  polar  fluid  by^^ 

pG^  =  gK-l  (12) 

Since  g^  increases  with  increasing  p  and  is  in  general  greater  than  one,  one  expects  that  L33  ( 
or  13^)  increases  with  increasirig  the  NLO  chromophore  density,  consistent  with  the  result  of 
Fig.  3.  In  ref.  7,  it  is  further  shown  how  the  dielectric  constant  of  the  combined 
chromophore/polymer  medium  affects  the  NLO  susceptibility,  due  to  the  orientational  pair 
correlation. 

The  orientational  pair  correlation  not  only  affects  the  NLO  susceptibility,  it  also  affects 
the  relaxation  behavior  of  the  NLO  susceptibility.  The  relaxation  behavior  is  directly  associated 
with  the  temporal  stability  of  the  nonlinear  optical  effect.  Shown  in  Fig,  4  are  the  representative 
temporal  characteristics  for  a  DAN  (10  wt  %)  doped  PMMA  film  at  68°C.  In  this  work  the 
sample  is  first  equilibrated  at  a  chosen  temperature  (68°C  for  that  of  Fig.  4)  then  the  poling  field 
is  applied.  After  the  SHG  signal  reaches  a  stable  maximum  value,  the  poling  field  is  then 
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switched  off  and  the  electrodes  are  shoned. 


As  one  notes,  the  SHG  signal  decay  behavior  cannot  be  fit  to  a  single  exponential.  The 
signal  first  assumes  an  initial  rapid  decay  and  is  then  followed  by  a  long  decay  to  the 
background.  As  shown  previously,  either  for  the  guest/host  or  for  the  chromophore 

functionalized  NLO  polymer,'^  the  SHG  signal  can  be  satisfactorily  described  by  a  fit  of  nonlinear 
optical  susceptibility  to  a  biexponential  form  as 

d,,  .  ixg  =  ae-'''  *  be*-  (13) 

where  a  +  b  is  equal  to  the  initial  value  of  the  nonlinear  optical  coefficient. 

Torkelson  et  al.  speculated  that  the  biexponential  decay  is  due  to  the  free  volume 
distribution  in  the  glassy  polymer,  consisting  of  liquid-like  and  solid-like  regions.  The  initial 
rapid  decay  is  associated  with  the  free  volume  distribution  in  the  liquid-like  region  in  which 
reorientation  of  the  chromophore  dipole  moment  is  rapid,  even  though  the  polymer  is  in  the 
glassy  state.  In  the  guest/host  system  we  found  that  the  initial  amplitude  decrease  depends  on 
the  strength  of  the  poling  field.  This  suggests  that  the  initial  decrease  is  probably  associated  with 
the  third  order  effect  similar  to  that  seen  in  EFISH  (electric  field  induced  second  harmonic 
generation).  The  time  constant,  tj,  characterizing  the  initial  decay  represents  an  effective 
relaxation  time,  which  also  convolutes  the  time  of  switching  off  the  poling  field.  However,  due 
to  the  slow  switching-off  electronics  used  in  the  present  experiment,  the  time  constant  rj  does 
not  reflect  the  true  third  order  effect  suggested  above.  Fast  electronics,  as  well  as  mode-locked 
laser  pulses  in  the  pico  or  subpico-second  duration,  are  needed  to  clarify  the  nature  of  the  fast 
initial  decay. 

Nevertheless,  the  second  slow  decay,  does  not  depend  on  the  poling  field,  but  strongly 
depends  on  temperature  as  well  as  the  chromophore  concentration.  The  slow  decay,  with  the  time 


Page  10 


constant  characterized  by  t,,  is  associated  with  the  orientational  relaxation  of  the  chromophores 
in  the  polymer  matrix.  As  the  temperature  increases  the  relaxation  decreases  the  polar  order  and 
hciicc  ;hc  .SiKj  stability.  We  shall  report  the  temperature  dependence  of  elsewhere  and 
consider  here  only  its  concentration  dependence. 

The  relaxation  time  extracted  from  fitting  the  nonlinear  optical  susceptibility  6^2  in 
accordance  with  Eq.  (13)  is  plotted  in  Fig.  5  as  a  function  of  NLO  chromophore  concentration. 
The  data  points  represent  values  obtained  for  the  film  at  68°C  and  poled  at  field  of  .89x10^ 
volt/cm.  Interestingly,  one  notes  that  increases  with  increasing  NLO  chromophore 
concentration,  suggesting  that  the  increase  in  the  loading  of  the  NLO  chromophore  density  slows 
down  the  SHG  signal  relaxation  and  enhances  the  stability. 

In  general,  doping  the  polymer  with  small  molecules  plasticizes  the  polymer  and  lowers 
its  glass  transition  temperature  Tg.  For  example,  doping  5  wt  %  of  DAN  in  amorphous  PMMA 
depresses  the  Tg  from  107  to  87°C.  Lowering  the  Tg  increases  the  mobility  of  the  NLO 
chromophore,  and  thus  we  expect  t2  to  decrease  with  increasing  chromophore  concentration. 
This  is  contrary  to  the  result  shown  in  Fig.  5,  which  indicates  that  T2  steadily  increases  with 
increasing  chromophore  concentration.  The  increase  in  t2  with  increasing  p  is  due  to  the 
orientational  pair  correlation  factor.  It  can  be  shown  that  if  represents  the  reorientational 
relaxation  time  of  an  uncorrelated  NLO  chromophore,  then  t2  is  related  to  by'^ 

72  =  73  (1  +  p  G2i)/(1  +  p  F)  (14) 

where  p  G^  represents  the  orientational  pair  correlation  factor  introduced  previously,  and  p  F 
represents  the  dynamic  orientational  pair  correlation  (DOPC)  factor.  This  situation  is  quite 
similar  to  light  scattering  in  which  the  recipiocal  of  the  linewidth  of  the  depolarized  Rayleigh 
spectrum  of  a  system  of  optically  anisotropic  molecules  is  proportional  to  the  orientational 
relaxation  time  7^.'^  The  orientational  relaxation  time  7j^  is  affected  by  the  static  and  dynamic 
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orientational  pair  correlation  in  a  manner  described  by  Eq.  (14).  In  the  present  situation,  we  are 
concerned  with  the  orientational  relaxation  associated  with  the  polar  order  L33,  rather  than  the 
relaxation  of  quadrapolar  order  involved  in  light  scattering. 

However,  it  appears  that  the  concentiation  dependence  of  tj  involved  in  SHG  is 
considerably  more  pronounced  than  that  observed  in  light  scattering.  Comparing  the 
concentration  dependence  of  73  with  that  of  L33,  one  notes  that  over  the  same  concentration 
range,  T2  increases  by  a  factor  of  more  than  4  (from  240  sec.  at  p  =  lxl0^‘^/cm^  to  1025  sec. 
at  p  =  3.7  x  10^®  (cm^))  in  contrast  to  a  factor  of  about  1.4.  The  larger  increase  in  73  indicates 
that  the  (DOPC)  factor,  as  given  in  Eq.  (14),  is  probably  negative,  which  decreases  the 
denominator  as  p  increases,  despite  the  factor  that  in  light  scattering  the  DOPC  factor  is  known 
to  be  small  and  often  negligible.  Clearly,  in  order  to  improve  our  understanding  of  the  effect 
of  orientational  pair  correlation  on  the  relaxation  behavior  of  the  SHG  signal,  additional  work 
is  needed. 
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Figure  Captions 


Fig.  1  The  boxcar  SHG  signal  output  signals  for  a  DNA/PMMA  film  at  84°C,  poled  with 
various  levels  of  the  electric  field  strength.  1;  ,61x10^  volt/cm;  2:  .77x10®  volt/cm;  3: 
.89x10®  volt/cm;  4:  1.04x10®  volt/cm;  5:  1.19x10®  volt/cm;  6:  1.31x10®  volt/cm 

Fig.  2  The  SHG  signal  (arbitrary  unit)  of  a  DNA/PMMA  film  (DAN  10  wt  %)  plotted  versus 
Ep^.  The  linearity  indicates  that  the  sample  is  poled  under  the  weak  poling  condition. 

Fig.  3  The  polar  orientational  order  parameter  L33  plotted  versus  number  density  (p)  of  DAN 
in  the  DAN/PMMA  film. 

Fig.  4  The  temporal  characteristic  of  the  SHG  signal  from  a  15  wt%  DAN/PMMA  film  poled 
with  an  electric  field  at  .89x10®  volt/cm. 

Fig.  5  The  relaxation  time  73  plotted  versus  number  density  p  of  DAN  in  the  DAN/PMMA  film. 
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